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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM NO. 1153 

CALCULATION OF THE "PRESSURE DISTRIBUTION ON BODIES 
OF REVOLUTION IN THE SUBSONIC FLOW OF A GAS 
PART I - AXIALLY SYMMETRICAL FLOW* 

By Herbert Bilharz and Ernst Holder 

I. INTRODUCTION 


The present report concerns' a method of computing the velocity 
and pressure, distributions on bodies of revolution in axially symmetrical 
flow in the subsonic range . 

The differential equation for the velocity potential $ of a 
compressible fluid motion is linearized in the conventional manner, 
and then nut in the form M = 0 by affine transformation. The 
quantity 0 represents the velocity potential of a fictitious incom- 
pressible flow, for which a constant superposition of ^.sources by 
sections- is secured by a method patterned after von Harman (reference 1) 

cS> 

which must comply with the boundary condition — » 0 at the origi- 

dn 

nally specified contour. This requirement yields for the "pseudo- 
stream function" a differential equation which must be fulfilled 
for as many points on the contour as source lengths are assumed. In 
this manner, the problem of defining the still unknown source intensi- 
ties is reduced to the solution of an inhomogeneous equation system. 

The pressure distribution is then determined with the aid of Bernoulli's 
equation and the adiabatic equation of state . 

Lastly, the pressure distributions in compressible and incom- 
pressible medium are compared on a model problem. 


"Zur Berechnung der Druckverteilung an Rotation skdrpern in der 
Unterschallstronung eines Gase3 Tell I: AchseTiRyrarr.eurischa Strbmung, " 

Zentrale fCtr wissenschaftliches Berichtsvesen tiler LuftfU^'tforschung 
(ZWB) Berlin- Adlershof, Forschungsbericht Nr. H69/I, Braunschweig, 

J an . 15 , 19^0 . 
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IT. STATIONARY POTENTIAL FLOW 


The discussion is restricted to a frictionless compressible fluid 
and it is assumed that the density of the medium is solely dependent 
upon the pressure p, while the functional relation "between pressure 
and density is to "be monotonic and continously differentiable. 

The steady motion follows then, for external conservative forces, 
the law of Lagrange analogous to the frictionless homogeneous incom- 
pressible fluid: If the frictionless compressible fluid moves 

irrotationally at time interval t Q , then it does so in every subsequent 
time interval t >t Q . Since the state of rest is a special case of 
the irrotational motion, the flow of a compressible fluid (nonstationary 
at first) starting from rest is free from rotation and certainly remains 
so as long as the flow velocity at no point becomes greater than the 
velocity of sound of the medium. An Irrotational motion, on the other 
hand, can always be represented by a velocity potential $, so that 
the velocity vector v 3 becomes the gradient of this potential. 


v = grad ^ 


( 1 ) 


with x, y, z denoting the rectangular Cartesian coordinates of 
the space (?) in which the fluid moves, and u, v, w the components 
of y in the directions of the axes of coordinates. Equation (l), 
expressed by components, reads 


t( 


3$ _ _ 
ax' " dy' 


hz 


( 2 ) 


On the assumption 


P * f(p) (3) 

the equation of motion of the compressible medium can be formally 
transformed into that of the incompressible medium by replacing in the 

latter, — by the ores sure function 

P 

p<p > = r§ <m 

5o 

is used throughout the test, and figures of this report in 
place of /»o, which was used in the original German report. 
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In the absence of external forces and in stationary flow, as invariably 
assumed in the following, Bernoulli’s equation reads 



+ P(p) 5 const. 


(5) 


Lastly, the equation of continuity is needed: 


idiv (pv) 


div v + 


I v o grad p = 
P 


0 


(6) 


The flow is then defined by (1), (5)# (6), and the subsequently 
discussed boundary condition (21)* 

Mow 


1 1 dp dp 

p s^ ad p = p 5p s rad p = 5 ^ s rad p 


hence by (6) in conjunction with (5) 


dp v2 

div v - v o grad 


= 0 


(7) 


Owing to 


v2 

v o grad “ 


7 o ^ X. ° 1L 




du 


+ TfC 


6y 


o dw 

ST 



+ 




equation (7) becomes, with observation of irrotationality, in 
coordinates. 
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with 



c 



the local velecity of sound* 

If the adiabatic equation of state 

p = p ”(; €? (?) 


is used as basis, equation ( 3 ) is satisfied; 
the pressure at any point of the medium and 
(undisturbed) flow at infinity; p and 

cp 

densities, K - r — is the ratio of specific 
c v 

mately IJ 4 I for air). By (5) and (O) the local sonic velocity 
point of the medium with the existing velocity is 


quantity p represents 
Pen the pressure in the 
are the corresponding 

heat of the gas (arprexi- 

at any 


and 


.2 - 




do) 



( 10 *) 


Therefore, o has a maximum value for y_ - 0, that is, at the stagnation 
point, where 


c max^ = 



JLl 

2 


lXaA 

oJJ 
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The critical velocity also follows from (10*) when c * |vf and 
resolved with respect to v: 


s crit 2 


c min 2 = 


2 

K + 1 





for 


isA = 0.8 and k = 1.1+08 
c oo 


for instance 


cmax = 1 * 0633^00 an ^ ciain = 0*9^9^20-0 


Since cmax is assumed at the contour while elsewhere in the 
medium 


Cmin < o < Cmax 


o may be put equal to in first approximation, by reason of 


22SJLZ-2E&. , 0.O887 
c 



0 , 201 ). 



0.0973 


narticular ly as c — with increasing distance from the contour 
as a result of the continuity of the flow process. 


III. LINEARIZATION 


Prom (8) it is seen that in all cases where all velocity 
components are negligibly small compared to the sonic velocity 
c min» equation (g) leads to equation 
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If the resulting velocity is no longer sufficiently small but 
still below Cmin» it is to be expected that the quantitative relations 
in the compressible medium are altered, but the character of the flow 
remains similar to' the incompressible medium. It is always assumed 
in the following that v2 < cmin^ in "the entire range of flow. 

Fith (£), equation (8) is written in the form 



and the values at infinity substituted in the coefficients of the 
second derivatives of 0 for u, v, w and c» Hence there follows 
by (12) a linear partial differential equation of the second order of 
the elliptical tyre with constant coefficients. 

If the x-axis of the coordinate system y points in direction 
of the velocity at infinity and 


v I u 

C < \ = r — < i 

°oo c oo 

denotes the ratio of undisturbed stream velocity to sonic velocity at 
infinity, the substitution 


T: u : u B , v = 0, w = 0 ; c = c M 


gives an appropriate linearization for (12). 

For the velocity potential 5 the equation reads 


(1 - \ 2 ) 


6x2 6y 2 6 z2 


( 13 ) 


It should be borne in mind that the linearization loses its 
validity near the stagnation points. 
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IV. PSEUDO -STREAM FUNCTION IN (7) 


Let 

-1 

A = + (1 - X 2 ) 2 
By affine transformation 

A: x = x, y = yA’ 1 , z ** zA" 1 

the system changes to 7. The equation 

1® (x, y, z) = $'(x, y, z) 

is to imply that to the points P(x, y, z) and "Pfx, y , z) coordinated 
hy A equal potential values are to correspond (transplanting of 
potential # from space (7) in space (7))* 

Owing to 

dx 2 Sx 2> by 2 Sy^ J 5z^ 

and (13), the potential 0 therefore satisfies the Laplace equation 




3 % 5 % 


Sz^ 


m) 


The transformation A therefore yields the rule: The flow 

around a contour C 2 in the compressible medium of a space (7) can 
he computed in linearized form in sin incompressible medium of a space (7) 
by reason of the fact that in (7) the contour 

C = AC 

affinely conjugated to C is used as basis. 


^2 is used throughout the text and figures of this report in 
place of the symbol ■cC' used in the orginal German report. 
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The following is restricted to axially symmetrical bodies and 
axially symmetrical flow; the x-axis is the axis of rotation. 

Putting , . , 


<? (3c, y) = <±> (x, y, 0) 


(15) 


we get, since the local function on each circle orthogonal to 
the x - axis has a constant value. 


Cp(x, r) = £> (5, y, z) * !?(*, r, 0) 


with 

r Z + (y 2 + z 2 ) 2 

and after some conversions 

r + = o, (? > 0) (16) 

6x2 6f2 

as the differential equation for the function <P • 

For this potential 9 there exists a function ^ , continuous with 
its first and second oartial derivatives and well defined up to an 
additive constant (not the Stokes stream function) which satisfies 
the equation 


6'j/ _ - 6cp 6 'ir _ - 6cp 

ox 6? 6f 6x 


The integrability condition for (17) reads 

6 / _6?N 6 / 6cp\ 

dr f dr/ = d2 V ) 


(17) 


( 18 > 
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or computed 



therefore (18) is fulfilled because of (16). By (17) then 


(x. 



By ( x 0 , r 0 ) is meant any point in the definition range of the function 
Cj)(x, r), that is, any roint of regularity of $(x, y, z) located on 
the olane z s 0. 

By (15) and (17) the velocity component of y. in axial and radial 
direction in a meridian olane (such as z * 0, for instance) is 



6$ _ dcp 
6* 6x 


]_Oj/ 
f dr 


(? > 0) 

- _ * _ 6$ _ 6cp _ i 6J 
dr 6? r dx 


(19) 


Introducing the terms resulting from (19) into the relation 

d / dcp~\ _ 6 /dcp\ 
dxVdr/ "dr\d5?/ 


finally gives the pseudo-stream function \|/ in ( y) : 
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V. BOTHDARY condition FOR 


Since C in ( y) is streamline, the boundary condition 


&£ - 0 

dn 


on C must be fulfilled for 


The derivation of the condition resulting *Yom (21) for 
in ( y) gives 


d$ d$ , x d$ , 

— - ; — cos (n, x) + — cos (n, r ) - 0 
dn dx br 


by (21)* 

Since 


b C' _ b Q ( ^)\[? 

Ox dx r dr 

6 $ i 6 $ 1 b\ \f 

or A d? Ar dx 


according to A and (1?) and since 


dr = dr_ dr. 
dx dr dx 


A 7“ - - ctg (n, x) 


we get 


Ad^dr 1_ d\£ _ 

r d? dx A r dx 


( 21 ) 


in * 
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hence 


A 2 ^ = 0 

6r dx 6x 


( 22 ) 


or 


A&2L+ o 

dr dr. dx 


— - « the direction of the tangent (as in the olane 2 = 0, for 
instance) to the contour. 

In corresrondence with the flow velocity (u^, 0) at infinity, 

^ r |uoof 2 + by (17). 

VI. CONSTANT AXIAL StI'-EHPOSTTIO BY SECTIONS 


After introducing three-dimensional T 'olar coordinates p, x 
in (y), that is, 

p length of vector 

$ angle between vector and axis x 

angle of meridian .plane with plane a " 0 

the velocity components are 


Up 


60 * c_ 
— , u^ 
op 


i 

p dT 


Ux 


0 


(25) 


The velocity potential for a simple source of yield Q in the 
origin is 
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On the "basis of the spherical symmetry the stream function 
is analogous to ( 19 ) 


60* _ 

l 

dilr* 

6 p 

P 2 sin § 

sO 

i 60* _ 

1 

_ 6\j/* 

p 6? " 

p sin 

% 6p 


hence 




_Q_ 

ilTT 


( 1 + COS §•) 


The stream function \j/ and the velocity component is now computed 
for a source of length a and constant yield q per unit length. 

The contribution of an element d£ to \j/ 0 in point P amounts to 


d\Jr Q z - -3. (l + cos &) dc, 
idTr 


hence \|/ Q is for the total length a 


r> a 


= " / (1 + cos $) d£ : 


~ (a + P* - P") 
4tf 


(21+) 


the latter because 


d^ cos ^ - dp 
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With 

Q : qa 

as the -tjotal yield of a, (21).) becomes 



(19) the axial and radial velocity components are 



(25) 


( 26 ) 


Because 


we get 


p s + (5^ + ?^) 


dp 

dr 


— = sin & 
P 


~E - * r cos 
dx P 



rh 


NACA TM Ho . 1153 


Hi 


hence finally by (26) 





O' 



Q 

i4.Tr a r 


( sin -V 1 


sin & 1 ) 


(cos 

i4.tr a r 


cos $.*) 



(27) 


Assuming consecutive source distances of constant length a and 
ield Qi (i = 1, 2, • • n ) on the x axis_ 1), and denoting 
he distances of a particular streamline point P from the ena points 
of the i-th source distance with p* * s and p\** =? pi + 1 
(i = 1, • • . . , n) we get by ( 25 ) 


n n 



\ 


( 28 ) 


This ^stream function 1 * is then sunernosed on that of the incompressible 
parallel flow (approach in positive 5 direction) 



so that the total flow is 


\lr - 




i=l 


( 29 ) 
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While the curves \{/ r constant in the_incompressible medium 
in (y) represent streamlines, the curve \J/ z 0 consists of especially 
the x- axis and 5, and \j/ z 0 must be applied for exactly as many 
points of 5 as there are constant source distances of length a with 
unknown intensity Qj- which give Qi by solution of a linear 
nonhomogeneous equation system* A somewhat different method must be 
applied in the compressible medium* 


VII. DETERMINATION OF SOURCE INTENSITIES 


The dimensionless quantities 



(30) 


are introduced as unknowns; = Pi k anc * Pj k" Indicate the 

length of the vectors reaching from the end noints of the i-th 
source distance to the k-th boundary point on C. By reason of 
the uninterrupted sequence of the source distances 


p i,k" = h 


+ 1, k 


( k, - 1 , 2, • 


n) 


The radius at the k-th contour point is indicated with r^; the boundary 
points on £ are to be so assumed that ?k (k = 1, * . . , n)is median 
perpendicular of the k-th source distance. 


Proceeding from (29) the linear equation system in the unknowns z± 
is setup. By it, in conjunction with, (30) 




n 



izl 


+ 



2 


a 


(5D 
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and by (27) 



(32) 


The boundary condition (22) for \j/ must be fulfilled from (31) • 
The final result with due allowance for (32) is 


_n 


sin & i+1 ) - 


1-1 



1-1 


- cos ^i+i) 


(33) 


This equation must be satisfied for n points on C, which have the 
distances (k = 1, 2 , • . . n) from the x axis: (33) is therefore 
written in the form 



c ik z 5 


a \dx/' k 


(k = 1, 2,, . n) 


(3U) 


r k and 
For the 


/dr \ 

Itajt ro ‘ er 

coefficients 


to the 
cik 


originally given contour 


S. i n 


(y). 
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°ik s (sin *ik ‘ sin ^i + l.k) 


+ 5 * (cos & i,k - cos ^i+l,k) 


( 35 ) 


Transforming (3l+) so that the new coefficients c*ik have the property 


s*kk *1 (k s 1, 2, . . «,n) 


(36) 


we get the equation system 


n 


^ c*ik zi 

/ 

i=l 


= | ■'(—) = (k = 

Vdx/k cos & kk 


■ 1. . . *»n) (37) 


with 


* - 1 
c*ik s c ik 


cos & kk 


where |cikl is small, when j i - kj is great. 
Putting 


ak = arc tan 


< 


dr^ 

•dx/k 


finally gives 


, _ * 

G ik 


1 . i i+l,k + ^i,k i . J ^i+l,k " ^i,k (, 7D , 

— slr \ - a k ) sinl 7(38) 

cos cos ^ ^ / 
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The formulas ( 36 ) and ( 38 ) prove themselves suitable for the numerical 
calculation by iteration. 


VIII. DETERMINATION OF VELOCITY DISTRIBUTION 


By (17) 

6\|/ _ 6cp 60 6\(/ _ _ 6cp _ . 60 

“ — H m, J* . mi. .... ■ . m v« — “ JT* — 

ox dr dr d r dx dx 


hence the axial and radial velocity comnonents on C: 


^ 6 5 i6)l ; w 5$ 

u x - r- = r , u r = — - = 

Ox r Or Or 


r dx 


(39) 


and the square of the magnitude of velocity 


I* - V * ^ * A**;*) . i(f f ♦ (f) 2 


(1+0) 


After computing the from (37)* and ^ p are known by (32), 

hence by (39) and (lj-0) the velocity component . 


IX. PRESSURE DISTRIBUTION 


By (10*) with the aid of (5) and (9) 


z/Jiy* 

vv 


W K-1 .2 A S 2 \ 

1 ~ l 1 - & 


\ 


2 


K 

in 


(i+i) ■ 
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and from it by expansion in powers of X^ (terms of power higher than 
the first in being disregarded); (See reference 2.) 


P_ 

£ 

2 




m 


while v 2 - jj2 

“CO 


and v_2 from (1+0) must be inserted in (41) and ( i+2) • 


Since the line of reasoning resulting in (42) holds only for 
small X, a different method fcr the pressure distribution is indicated, 
which holds for all X with 


0 < X s - < max 

— c oo 



and which gives the Mach number at any point of the body and from it 
the local velocity of sound c. 

The quantity w is introduced through the equation 



-2JL.BB.* v 2 

K - 1 Pco 


by (43) 


w proves itself as maximum velocity; and (10) becomes then 


5 2 = JSfi (1 - I?) 


(It3) 





c2 = 


c 2 

«2 


with 
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with c replacing c, hence putting jw J r 1 the square p. of the 
Mach number is by (1+3) 


_ _ 2 TSf 2 

^ c 2 K - 1 1 - w2 


m 


and therefore 


jw| 


+ 


(k - 1)|J. 

2 + ( K - ] )p 


1 

2 


(U5) 


From (ii) and (9) follows 


dp r 


dp 

P 


Co p K " 2 d p 


and from it by integration 


P = Cl p K_1 : Cl 


(U6) 


KPo 

C 0 - ? and Cl are constant. 

Po 

By (1+6) 


,2 = 


dp _ dP 
dp d log p 


= (K - 1) P 


(1+7) 


-V 2 + P 

2~ 



m 


On the other hand 


by (5) 
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••2 

hence, heoause of i : 1 

P = \ (1 - w 2 ) 

Therefore by (1+7 ) in conjunction with (1+3 ) 

JL-L-I (l - s 2 ) : c 2 : ( K - 1) P = (k - 1) C x p 

2 

hence 

Cx P K " 1 = \ (1 - w 2 ) 


with p 0 as the density in the stagnation point the integration 
cons bant Ci becomes 



1 - K 


hence by (U9) finally 



1 - 


w 2 


If p 0 is the pressure in the stagnation point 
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by (50), whence by (I4I4) 


Po _ 
P 




K 

K-l 


therefore 


P - 


K-l 


K-l 


Po 

VPv 


V- 


(51) 


Because c2 



we get by (51) 


W' 


2 = q 2 


P - 


2k p 
e - 1 P 


K-l 

fV\ K 

Uv 


- 1 


hence 



(52) 


Let 



(53) 
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that is finally 



the value from to be inserted for n* 


X. EXAMPLE 


The foregoing considerations are suggestive of exploring for 
the first those bodies of revolutions the meridian curve C of 
which is an analytical function of x in the (x, r) plane: 

C r r = r(x) 


The example proceeds from the series of ourves 


(n + 1)“ A 

r = d £3 x a ( 1 - x) 

with 


d>0, 0<n <1, 0<x<l 


The factor d gives a measure for the bhiokness of the body. 

Tf 0 <n< 1 the profiles in the (x, r) plane have a round forebody 
and a pointed afterbody. For n z 1, C is symmetrical to x = 0.5* 
in which case the forebody and the afterbody are pointed. 

Generally 


r} max = d for x = 

n + 1 
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The values chosen for the parameters are 

d as n = \ = 0.8 (hence , Um = 0.64)) A = + (1 - X^) 2 - 2 

Since r = A"- 1 - r = 0.6 r, the meridian curve 0 of the affine trans- 
formed body lies within C. The meridian sections g and c are 
shown in figure 3. 

The constant length a of the source distance was put at O.i. 
After computing the coefficients ( 38) , the system (37) is solved 
hy iteration. Thereby it was found that the values Zj_ for the 
source intensities converge quicker than hy the von Karman method for 
the initial hody in incompressihle How. Both distributions of 
vortex intensity for the same body of revolution in compressible and 
incompressible flow for the same flow velocity at infinity 

U = 266 .4m/ sec" 1 , u<» = 0.64 are included in figure 3* 

Figure 4 shows the velocity in terms of U) figure 5 the Mach 
number for compressible flow. The Mach number follows from (4)4) as 

function of the velocity |wj = |w M j, with jw» j = 0.3363 
according to figure 2 and equation (45) . 

Figure 6 represents the pressure distribution by ( 54) on the body 
along with the pressure distribution of incompressible flow for 
comparison. 

The numerical values are given in table I. 


XI. SUMMARY 


The present report gives a method for determination of the 
velocity and pressure distributions as they appear on bodies of 
revolution in axially- symmetrical flow in the subsonic range of a 
gas with adiabatic equation of state . 

The differential equation for the velocity’- potential $ of 
the gas motion is first linearized in such a manner that in the 
coefficients of the second partial derivatives of $ with respect 
to the coordinates for the velocity components u, v, w and the 
local sonic velocity c the corresponding values from infinity 
U = tU, 0,0j Co, are inserted. There results for <3? a linear 
partial differential equation of the second order of elliptic type 
with constant coefficients. By an affine transformation A of the 
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variables it can. be converted to the Laplace form £0 = 0. The 
performed linearization is valid only outside of a certain 
neighborhood of the stagnation points. - 

By the transformation A the body of revolution with the 
meridian curve C:r = r(x) is transformed into a similar body 
normal to the axis of rotation with the meridian curve « r( x) . 

$ as solution of the equation M - 0 represents the velocity 
potential of an incompressible fictitious flow in a space (y ) j it 
is to be noted that £T is not a streamline. On the contrary,- the 
equipotehtial lines 0 = Constant intersect the curve _C not 
orthogonally. Thus $ is to be determined as a^regular potential 
f motion which satisfies in the outer range of C the equation 
M - 0, _at infinity has the form + reg. potential function, 

and on C fulfills the condition 


1 SO dr- 

2 c)x dx“ 



0 


This last requirement results from conversion of the boundary 
condition SO/Sn = 0 to Q for the velocity potential 0 . 

As in the incompressible medium there exists for 0 a function 
"pseudo- stream function" which satisfies the conditions 

-r M • M M 

dx Sr’ ^r Sx 


in the outer range of C[ satisfies the equation 




at infinity has the form \jT e= 1 Ur 2 + regular function, and on € 
satisfies the equation: 2 
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with 



-I - 

M-oo) 2 


Under the assumption that the function $ can "be analytically con-, 
tinued into the interior of the body to the symmetry axis, $ can 
be produced by a constant superposition of sources on the axis. 

Ton Karman's method mentioned at the beginning is based on this 
assumption: the continuous source distribution on. the axis is 

replaced by a constant superposition in sections. The number of ■ 
source lengths is selected according to the desired degree of • 1 
accuracy. Corresponding to this number n one chooses on the 
contour £ n varying "control points j" most serviceably, on the 
center perpendiculars of the source lengths . The solution of the 
integral equation for the continuous source distribution is reduced 
to a linear inhomogeneous equation system of the nth order, the 
unknown source intensities determined so that they satisfy the 
equation (55) in the control points. 

The solution of this equation system by iteration converges 
more rapidly for a body in compressible medium than for the same 
body in a fluid of constant volume, the reason is on the one hand 
the form of the coefficients which is more symmetrical than von Karman' 
on the other hand the n-eater slenderness of the auxiliary body (c) 
which results by the transformation A. 

The pressure distribution is determined from Bernoulli ' s 
equation and the equation of state. It is advantageous to introduce 
the limiting velocity w (outflow velocity of the gas into the 
vacuum) . Thus it is possible to determine the Mach number at every 
point of the contour and to compare pressure distributions for 
various Mach numbers. 


As example the pressure distributions for a certain body of 
revolution were calculated in the compressible medium (|i ra = 0.64) 
and in the medium of constant volume. The results are compiled 
in table I. 


Translated by J. Tanier 
National Advisory Committee 
for Aeronautics 
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TABLE I 


i 

X 

;r| 

d | r | 

z, compressible 
1 

Zj incompressible 

L r l compressible 
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